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Interest in epitaxial ferroelectric nanoislands 
was growing rapidly in recent years driven by 
their potential for devices, especially ultradense 
memories [fl |2l [Sj I2- Recent advances in 
the "bottom-up" (self-assembly) nanometer scale 
techniques have opened up the opportunities 
of fabricating high-quality epitaxial ferroelectric 
nanoislands with extremely small thickness and 
lateral size on the order of 1 nm and 20 nm, re- 
spectively [41 m m m l]. On the other hand, 
recent emergence of powerful probes, such as 
piezoresponse force microscopy (PFM), has en- 
abled imaging of a local domain structure with 
sub-10 nm resolution [3, 9, lOj. In spite of those 
developments, a clear understanding of the polar- 
ization patterns in epitaxial ferroelectric nanois- 
lands is lacking, and some important character- 
istics, like a critical lateral size for ferroelec- 
tricity, are not yet established. Here, we per- 
form ab-initio studies of non-electroded epitax- 
ial Pb(Zro.5Tio5)03 (PZT) and BaTiOg (BTO) 
nanoislands and show the existence of novel po- 
larization patterns driven by the misfit strains 
and/or anisotropy energy. The results allow in- 
terpretation of the data and design of the ferro- 
electric nanostructures with tailored response to 
external field. 

Geometrically, nanoislands or flat nanoparticles rep- 
resent a class of systems that bridge the gap between 
OD nanodots and 2D ultra-thin films. Compared with 
thin films, they have free-standing side walls that tend 
to suppress the formation of uniform in-plane polariza- 
tion because of appearing depolarizing field, like in fer- 
roelectric [lU[T2] or ferromagnetic particles [13 . On the 
other hand, relative to the (confined in all three dimen- 
sions) nanodots, they have large aspect ratio and likely 
to behave similarly to thin films when the polarization 
is out-of-plane. These observations lead one to expect 
that ferroelectric nanoislands should exhibit some kind 
of a "particle-to-thin film" crossover behavior and related 
novel effects depending on the aspect ratio and the type 
of bulk polarization ordering (which is different in PZT 
compared to BTO crystals). 

Free-standing ferroelectric nanodots usually lose sta- 
bility with respect to a vortex ground state, in which the 
polarization curls around some vortex core(s). Such an 
ordering can be characterized, in simplest single vortex 
cases, by a toroidal moment or a moment of polarization 
[m [T5J [ini HZ] . This is analogous, to some extent, to the 



magnetic vortices in magnetic particles |13| . but there 
are also qualitative differences due to electrostatic inter- 
actions and inhomogeneity energy in ferroelectric parti- 
cles being larger than the magnetic analogues by many 
orders of magnitude. Note also that the symmetry of 
the bulk cubic crystals of BTO or PZT does not allow 
for a "bulk" toroidal moment |18j . such an ordering in 
the present case is entirely due to presence of bound- 
aries of the particles and the accompanying it depolar- 
izing/demagnetizing field. The moment of polarization 
appears continuously at the critical point and may be re- 
garded as an order parameter for such a transition in an 
island [Tij . 

We shall consider below the nanoislands of cubic per- 
ovskite materials. The epitaxial thin films of such mate- 
rials become uniaxial under compressive strain with an 
easy axis perpendicular to the film, and usually form a 
180° stripe (c-)domain structure [11 [2S IHl [H HI [Ml 
E51 E51 \T7\ E5] . One can, therefore, anticipate the flat fer- 
roelectric nanoislands with large aspect ratio to undergo 
an unusual vortex to stripe domains transformation with 
an increasing compressive misfit strain. We show below 
that perovskite nanoislands do indeed undergo such a 
transition, but only via one or two intermediate phases 
where in-plane curling polarization coexist with the 180° 
c-domains. Moreover, the domain phase does not neces- 
sarily render the well known 180° stripe domain struc- 
ture, but also can take the form of a 180° tweed texture, 
depending on the anisotropy energy. We also demon- 
strate that stripe and tweed textures evolve differently 
in applied external electric field. 

The simulated nanoparticles have circular and/or rect- 
angular shapes and the following dimensions: 19x6, 
39x 10 for the former and 39x39x 10 for the latter (all di- 
mensions are given in lattice parameters of the cubic bulk 
phase a, which is 4.00A and 3.95A for PZT and BTO, 
respectively). The z— axis is selected along the (shortest) 
pseudocubic [001] direction normal to the substrate. To 
simulate finite-temperature behavior of the particles, we 
use the effective Hamiltonian 

H{u,} = H''{u,} + H'{u,} - J2 Z*u, ■ Eo, (1) 

i 

whose primary variable is Ui, the local soft mode in the 
unit cell i [SSJ |3D]. The first part of the Hamiltonian, 
H'' is the same as in infinite bulk materials: its first- 
principles-derived parameters for BTO are given in [29], 
and for PZT- in [30]. The second part, , represents 
the corrections to the H'' associated with the presence 
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FIG. 1: Low-temperature polarization distributions on the central z-plane in PZT (a,b,c) and BTO (d,e,f) nanoparticles in 
free-standing states. (a,d)- 19x6 nanodisks, (b,e)- 39x10 nanodisks and (c,f)-39x39x 10 particles with square footprints. 
Shown in red, blue, green and black are the dipoles directed predominantly along [±100] and [0±10] in PZT and along [±1±10] 
in BTO systems. 



of free surface via "vacuum-local mode" and "vacuum- 
inhomogeneous strain" interactions El]. And, fi- 
nally, the third term describes the interaction of the lo- 
cal site dipoles = Z*Ui with an external field Eq [33] 
{Z* is the effective Born charge of the local mode M;). 
We should stress that the first-principles-derived Hamil- 
tonian is able to reproduce some nontrivial properties 
beyond the properties it was fit to: phase diagram, oc- 
currence of an exotic monoclinic phase in a small range 
of Ti concentrations in bulk PZT, stripe c-domains in 
thin films with compressive in-plane strains, etc. Typ- 
ically, 40,000 Monte Carlo sweeps are used to find the 
equilibrium dipole configuration at a fixed temperature. 
Open-circuit boundary conditions are assumed, so that 
no screening charges are taken into consideration. The 
effect of a substrate is imposed by fixing the homoge- 
neous in-plane strain in the whole island: e^x = £yy = e, 
exy = 0. 

Our calculations predict that at low temperatures and 
in free-standing states all the investigated nanoparticles 
transform into a vortex ground state characterized by a 
nonzero toroidal moment G = {2N)~^ J2i '"i ^ Pi> where 
N is the total number of unit cells. In disk-shaped PZT 
particles (Fig. la,b), the local dipoles rotate from cell to 
cell forming the classical vortex patterns. But in a rect- 
angular PZT 39x39x10 particle the polarization aligns 
with the square boundary as much as possible, hence the 
discontinuity lines (domain walls) form, and the resulting 



pattern can be viewed as consisting of four domains sep- 
arated by 90° domain walls (Fig. Ic). This is in a clear 
analogy with ordering in magnetic particles |13j . antic- 
ipated for ferroelectric particles back in 1950s [TTJ fT2] . 
The polarization in each domain point along one of the 
four pscudocubic [±100] and [0±10] directions, so that 
the domain walls are parallel to the (1±10) crystallo- 
graphic planes and pass through the center of the par- 
ticle. In BTO particles, regardless of their shape and 
size, the dipoles tend to be directed along the [±1±10] 
directions (Fig. ld,e,f). As a result, two (lOO)-type 90" 
domain walls are likely to be formed crossing each other 
at the geometric center of the particles. The walls are 
rather fuzzy in bigger particles, especially in a rectangu- 
lar 39x39x10 dot where the dipoles are frustrated near 
the lateral surfaces. 

The different behavior of PZT and BTO free-standing 
nanoparticles can be easily understood by looking at the 
anisotropy energy J2^ lA^lx^ly + "LuL ± u\u\^) in the 
effective Hamiltonian H (its part iJ*"). Suppose that the 
amplitude \u\ is fixed, then the anisotropy energy will 
approach its minimum either for u along the [100] direc- 
tion (7 > 0) or along the [111] direction (7 < 0). For 
PZT, the parameter 7 (associated with the dominating 
Ti-centered unit cells) is positive [30] - this is consistent 
with the fact that PZT in bulk form adopts a tetragonal 
structure with an "easy" direction of polarization along 
[100]. At the same time, in BTO 7 < [21], and this 
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system in bulk has a rhombohedral equilibrium structure 
with easy direction [111] and the equivalent bulk diago- 
nals. Let us compare now, for example, two 39x39x10 
PZT and BTO particles. In both cases, the effects of de- 
polarizing field force the dipoles to lie in the x, y— plane. 
In the case of the PZT particle, however, the materials 
anisotropy factor is in compliance with the shape: ma- 
jority of the dipoles are oriented along the easy directions 
[100] and [010] and at the same time they are parallel to 
the surface. The situation is quite different for the BTO 
particles, where none of the dipoles can be directed along 
the easy [111] body diagonal: in this case the majority 
of dipoles in the center region prefer to be oriented along 
[110] instead of the [111] direction. Near the side surfaces, 
however, the depolarizing field forces the dipoles to align 
with the square geometry, which increases the anisotropy 
energy term. The competition between these forces leads 
to the dipole frustration near the vertical walls and to 
the formation of four additional vortices near the corners 
(Fig. If) in agreement with the previous calculations for 
a 24x24x24 BTO nanodot 

The tensile strains stabilize the vortex states further, 
as seen from Fig. 2, increasing the ^-component of the 
toroidal moment. On the contrary, compressive strains 
enhance the out-of-plane polarization, decrease Gz and 
eventually lead, as analysis below shows, to a 180" multi- 
domain structure running through the entire thickness of 
the islands. It is remarkable that in the case of BTO 
particles the critical strain associated with this transition 
is very small, approximately zero, ~0 %. In the PZT 
case, the critical strain dramatically decreases from 1.2 
% to «0 % when one increases the diameter from 19 to 
39. The bigger PZT particle in a vortex state can tolerate 
the compressive strains less because its optimal in-plane 
averaged lattice parameter is noticeably larger than that 
corresponding to the bulk cubic phase (Fig. 2b). 

In the PZT nanoparticles, the 180° multi-domain 
phase is nothing but the well-known 180" stripe domain 
structure found previously in infinite PZT thin films 
[SgEl]. Here, like in the thin films, the "up" and "down" 
domain run along [100] and alternate along [010] direc- 
tion (Fig. 3a,b). The domain width, however, is not ho- 
mogeneous across the particles: it is larger in the interior 
region and becomes 1-2 unit cell narrower near the edge 
of the particles. Besides, the stripes become wider in 
bigger particles. In a circular 6x19 particle, for example, 
the width of the central stripe is 4 unit cells (16 A), while 
it is 6 unit cells (24A) in 10x39 and 10x39x39 particles. 
In both cases the stripes are narrower than that found in 
the PZT thin film modeled with a 40x40x5 supercell (8 
unit cells), despite the fact that our particles are thicker. 

BTO nanoparticles, in comparison with PZT coun- 
terparts, adopt not straight, but curved zig-zag-type 
stripes, running predominantly in [110] and [1-10] direc- 
tions (Fig. 3c, d). Though rarely, they intersect each other 
forming a 180" tweed domain texture, which is drastically 
different from that found in BTO thin films pUll^ . In 
the case of thin films, the stripes run only along one di- 
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FIG. 2: The z-component of the toroidal moment as a func- 
tion of strain in 19x6 (a) and 39x10 (b) nanoparticles. 
The strains are measured relative to the theoretical LDA- 
calculated lattice constants in bulk cubic structures. Note 
that in free standing states the particles are characterized 
by non-zero "residual" strains marked by vertical lines. The 
residual strain increases in passing from 19x6 to 39x10 sys- 
tems, especially in PZT. 



rection (either [110] or [1-10]) and the system loses its 
four-fold symmetry axis perpendicular to the film. In 
the nanoparticles, however, the axis of the 4-th order is 
(on average) preserved by allowing the forming stripes to 
propagate in both possible directions. The stripes width 
changes from place to place, being less than 4 lattice con- 
stants on average. This value can be compared with 4.3 
lattice periods found in BTO thin films [28]. 

As the present analysis shows, the transformation 
"vortex-to-180" domains" does not occur directly, but 
rather via one or even two intermediate structures. The 
simplest case is presented by a 19x6 PZT particle having 
only a single intermediate structure. This structure de- 
velops from the vortex core that turns out to be longitu- 
dinally polarized already in a free-standing state. While 
the core region (with the radius 3-4 lattice spacings) em- 
braces the polarization pointing "up" , beyond this radius 
the polarization tilts weakly out of plane oppositely to the 
core region, so that sum of all the dipoles add up to zero 
(Fig. 4). Under compressive strains, the vortex core ra- 
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FIG. 3: Low-temperature three-dimensional polarization 
patterns in free-standing (001) PZT (a,b) and BTO (c,d) 
nanoparticles under compressive strain of -2%. (a,c)- 39x10 
nanodisks, (b,d)- 39x39x10 particles with square footprints. 
The small rectangular inserts show the cross sections of a 
central vertical plane indicated by yellow lines. Red and 
blue arrows show local dipoles with positive and negative z- 
component, respectively. 



dius expands and becomes comparable to the lateral size 
of the particles. At this stage, the x, y— components of 
the local dipoles still form a vortex, while the out-of-plain 
component, z, breaks the system into coaxial oppositely 
polarized cylindrically symmetric domains-skyrmion-like 
structure [M]. As the compressive strain further in- 
creases, this structure abruptly transforms into a 180'^ 
stripe domain phase with zero moment G (Fig. 2a). 

In all other particles (both PZT and BTO), the 
structural transformation into 180" domain phase oc- 
cur in a two-step fashion. At first and relatively short 
stage, again, two symmetry conforming 180° coaxial do- 
mains develops at the expense of the toroidal moment 
(skyrmion structure). As the in-plane compression in- 
creases, this skyrmion structure transforms into 180° 
stripe (tweed) domains, as demostrated in Fig. 5 for a 
39x10 PZT nanoparticle. Remarkably, this happens on 
the background of nonzero Gz- And finally, the toroidal 
moment disappears at a critical compressive strain, and 
the system enters into a pure stripe or tweed domain 
phase. 

If the compressive strains are released, the system in a 
180" domain state relaxes back to a vortex dipole config- 
uration. Typical evolution of the toroidal moment and 
the average value of the z-component of the local mode 
1^2 1 — N~-^ Wi,z \ with the number of the Monte-Carlo 
sweeps is shown in Fig. 6. It is seen that the Gz compo- 
nent first linearly and then non- linearly increases, until 
saturation is reached, the other two components being 
practically null. The \uz\, which is simply proportional 



FIG. 4: Vertical polarization component in a 19x6 PZT disk 
at different misfit strains. While the values are radially av- 
eraged, they correspond to an instant moment of time and 
represent the fluctuating dipole structure. The lines repre- 
sent 4-th order polinomial fits to the instant polarizations. 




FIG. 5: Distribution of polarization in a 39x 10 PZT nanodisk 
for two close strains of 0.6 % (left) and 0.5 % (right). The 
latter, which is more compressive relative to the free stand- 
ing nanoparticle, leads to a state where 180" domains coexist 
with the curling polarization. The small rectangular inserts 
show the cross sections of a central vertical plane indicated 
by yellow lines. Red and blue arrows show local dipoles with 
positive and negative z-component, respectively. 

to the average magnitude of the out-of-plane polarization 
\Pz\, decreases in a similar way, until it becomes satu- 
rated at some non-zero value. The last fact is mainly 
connected with the oppositely polarized core ( "up" ) and 
periphery region ( "down" ) . 

Consider now how the stripe domain and tweed struc- 
tures evolve under the external electric field. The typi- 
cal behavior of PZT systems is presented by a cylindrical 
particle 10x39, whose evolution under different fields ap- 
plying along z is shown in Fig. 7. On can see that the 
dipoles with the direction opposite to the applied field 
start fiipping first near the particle edge. When the field 
is further increased, the "switching front" propagates to 
the center, shrinking the length of the stripes and (to a 
lesser extent) their width, still being opposite to the field. 
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FIG. 6: The evolution of the toroidal moment G and average 
magnitude of the z-component of the local mode in a 39x10 
PZT disk after "switching" off compressive strains. 

Ultimately, one reaches a state in which only one reduced 
stripe remains, near the particle center; this stripe has 
an eUiptical cross section and can be called a "bubble" 
(Fig. 7, S = 50 X 10^ V/cm). BTO nanoparticles be- 
have qualitatively similar to the PZT counterparts only 
in the initial stage of evolution: there the dipoles also 
start flipping in the field direction in the rim region caus- 
ing the longitudinal shrinking of the stripes. Now, how- 
ever, the stripe shrinking is accompanied by their break- 
ing into pieces (substripes) which, in turn, transform into 
nanobubbles. As a result, several relatively thin bubbles 
form near the center before the systems enters a mon- 
odomain state (Fig. 7, 40 x 10^ V/cm). Note that both 
in PZT and BTO, the out-of-plain polarization linearly 
increases with the electric field, similar to the case of 
ultra-thin films [?7l[^ . 

All the considered above intermediate phases (includ- 
ing skyrmion-like) are characterized by z-domains coex- 
isting with the in-plane curling polarization; they can be 
called the "domain-patterned vortex" structures. Inter- 
estingly, such structures have never been reported for 
ferroelectric small systems, though has been found in 
flat magnetic particles with low perpendicular anisotropy 
|35j . Analogously, the 180° tweed domain texture has 
never been discovered in ferroelectric nanostructures, 
though resembles pretty much the labyrinthine phase ob- 
served in ultrathin magnetic films |36j . One can think 
that this structure is metastable and eventually will 
evolve into an ordinary stripe phase. In our simulations, 
however, we did not observe any signs of such transitions 
even after a huge running time, more than 500 000 Monte 
Carlo sweeps. 

Our predictions can be helpful for understanding some 
puzzling and controversial experimental results 5, 6, 7, 
|B] . While the authors [51 17] did not observe any piezore- 



sponse of the PbTiOs islands with the diameters smaller 
than 20 nm, Ref. [6] reported observation of ferroelec- 
tricity in zirconate titanate nanoparticles with the lateral 
size as small as 9 nm. On the other hand, according to 
measurements [F], the relationship between the size of 
nanoislands and occurrence of ferroelectricity is rather 
random, raising the question why some particles are fer- 
roelectric, but the other (even with very similar sizes) 
are not. Based on our results, we can explain these un- 
certainties in the experimental results in the following 
way. As compared with continuous thin films, the exis- 
tence of free side walls in nanoislands leads to new and 
additional mechanisms of the strain relaxations |37j . in- 
cluding formation of partial or full misfit dislocations at 
nodes of the lateral surfaces and the fiat plane of the 
substrate. This implies that with the shrinking in lateral 
size, the maintaining of coherent mismatch lattice strains 
in a nanoisland becomes more and more difficult, and it 
is very likely that under the same experimental condi- 
tions the smaller particles will experience less external 
strains. Moreover, since the magnitude of strains is not 
well-controlled, even the particles with comparable sizes 
can be differently strained. But in this case, according to 
our calculations, they can show totally different piezore- 
sponses. Indeed, as the compressive strain is relaxed, an 
island can transform into an intermediate or even in a 
vortex state, which is unable to produce any (vertical or 
lateral) PFM signal at all. The particle in a vortex state 
should be viewed as a paraelectric one showing only linear 
and hysteresis-free piezoresponse as a function of dc bias 
voltage applied to the AFM tip. Similar arguments can 
be used in interpretation the results of Ref. [31] report- 
ing the absence of the local piezoresponse loop in a PZT 
nanoisland with a height of 10 nm [38 . Although the au- 
thors attributed this effect to the inhomogeneous elastic 
fields associated with misfit dislocations [38], it also can 
be explained by the formation of a curling polarization. 

In summary, we investigate the strain effects on (001)- 
epitaxial PZT and BTO nanoislands using first-principle- 
derived effective Hamiltonian approach. Our study leads 
to the following conclusions: (i) regardless of size, shape, 
and material, all the investigated particles transform into 
a vortex state with in-plane curling polarization; the lat- 
ter resembles classical vortex pattern in smaller/round 
particles, while in bigger/rectangular particles it can also 
be described as a domain structure with four closure 
domains and four 90*^ domain walls, (ii) under strong 
enough compressive strains the vortex state is no longer 
stable and yields to a 180° stripe domain phase in PZT 
and to a more symmetric, 180" tweed domain texture, 
in BTO nanoparticles, and (iii) these transitions proceed 
not directly, but via one or two unusual intermediate 
phases where domain structure coexist with circular vor- 
tex ordering. We further discovered that the electric-field 
induced evolution (switching) of multi-domain structures 
is essentially different for the two materials (PZT and 
BTO). 
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